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f~ s , Abstract 

A special class of symmetry reductions called nonclassical equivalence transformations is discussed in 
connection to a class of parameter identification problems represented by partial differential equations. 
These symmetry reductions relate the forward and inverse problems, reduce the dimension of the equa- 
tion, yield special types of solutions, and may be incorporated into the boundary conditions as well. As 
an example, we discuss the nonlinear stationary heat conduction equation and show that this approach 
permits the study of the model on new types of domains. Our MAPLE routine GENDEFNC which uses the 
package DESOLV (authors Carminati and Vu) has been updated for this propose and its output is the 
nonlinear partial differential equation system of the determining equations of the nonclassical equivalence 
transformations. 

Keywords: classical Lie symmetries, symmetry reductions, parameter identification problems 
AMS: 58J70, 70G65, 35R30, 35G30 

1 Introduction 

\Q , One of the fastest developing research fields in the last few years is the area of inverse problems. These 

problems arise frequently in engineering, mathematics, and physics. In particular, parameter identification 
problems deal with the identification of physical parameters from observations of the evolution of a system 
and especially arise when the physical laws governing the processes are known, but the information about 
the parameters occurring in equations is needed. In general, these are ill-posed problems, in the sense that 
they do not fulfill Hadamard's postulates for all admissible data: a solution exists, the solution is unique, and 
the solution depends continuously on the given data. Arbitrary small changes in data may lead to arbitrary 
large changes in the solution [TT] , 

Consider the following class of parameter identification problems modeled by partial differential equations 
(PDEs) of the form 

<3 ; F(x,u> (m) ,£; (n) ) = o, (i) 

where the unknown function E = E(x) is called parameter, and, respectively, the arbitrary function w = w(x) 
is called data, with x — {x\, . . . , x p ) G Vt C R p a given domain, where w^ m ' denotes the function w together 
with its partial derivatives up to order m. Suppose that the parameters and data are analytical functions. 
The PDE |T]) sometimes augmented with certain boundary conditions is called the inverse problem associated 
with a direct (or forward) problem. The direct problem is the same equation but the unknown function is 
the data (for which certain boundary conditions are imposed as well). 
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On the other hand, symmetry analysis theory has been widely used to study nonlinear PDEs. A re- 
markable number of mathematical physics models have been successfully analyzed from this point of view. 
There is a considerable body of literature on this topic (see, for example, [T], @]-[5], [S], PU], [T3]-[T5]. 
[TS]-[2n], E2|-[23 and references from there). Sophus Lie [T7] introduced the notion of continuous groups 
of transformations (known today as Lie groups of transformations) and, subsequently, a method for finding 
the symmetry group associated with a PDE. This (local) group of transformations acting on the space of 
the independent variables and the space of the dependent variables of the equation has the property that it 
leaves the set of all analytical solutions of the PDE invariant. Moreover, the form or the equation remains 
unchanged and group-invariant solutions can be found. Lie's method has been applied extensively to various 
mathematical models represented by PDEs. Consequently new methods for seeking explicit solutions that 
cannot be obtained by Lie's method have been proposed over the years. The nonclassical method intro- 
duced by Bluman and Cole : A\ is a technique for determining the nonclassical symmetries (or conditional 
symmetries) related to a PDE. These are transformations that act on the space of the independent vari- 
ables and the space of the dependent variables with the property that they leave only a subset of the set 
of all analytical solutions invariant. Knowledge of these classical and nonclassical symmetries allows one to 
reduce the dimension of the studied PDE and to determine solutions of special type. It may happen that 
certain group-invariant solutions cannot be found explicitly, but at least additional information about the 
studied model can be obtained in this case. Any classical symmetry is a nonclassical symmetry but not 
conversely. Ovsyannikov [5U] introduced the notion of equivalence transformations for a PDE depending 
on an arbitrary function and these transformations have been recently generalized by Meleshko in 18J. 
The (generalized) equivalence transformations are groups of transformations acting on the space of the in- 
dependent variables, the space of the dependent variables, and the space of the arbitrary functions with the 
property that they leave the differential structure of the equation unchanged. By dropping the auxiliary 
conditions that characterize the functional dependence, Torrisi and Tracina introduced the notion of weak 
equivalence transformations |24j . These transformations have been later generalized by Romano and Torrisi 
in [23]. Notice that all these methods do not take into account the boundary conditions attached to a PDE. 
Moreover, in the equivalence transformations approach the arbitrary functions may depend on the dependent 
variable only. In our case, for parameter identification problems of the form (fTJ), the data does not depend 
on the parameter - both data and parameter act on the space of the independent variables. 

The aim of this paper is to continue the study of the parameter identification problems ([T]) from the 
point of view of symmetry analysis theory. This is a second paper in which we analyze symmetry reductions 
suitable for (fl]) . Investigation of special Lie groups of transformations related to a model may allow us to 
reduce its dimension and further obtain information about its solutions. As discussed in [3J, finding the 
classical Lie symmetries related to the inverse problem (flj might be a very difficult task and the success 
of the method depends on the nonlinearity of the equation. Therefore, considering both the parameter 
and data as unknown functions in (fT]), the classical Lie method can be applied with less difficulties (in [3] 
a mathematical model arising in car windshield design was analyzed using this assumption). Throughout 
this paper, similar to [3J, we shall refer to these transformations as equivalence transformations, taking into 
account that w was initially an arbitrary function (these transformations are different than the equivalence 
transformations introduced by Ovsyannikov, but in the literature they are also referred to as equivalence 
transformations). Furthermore, once the equivalence transformations related to ([T]) have been discussed, 
it is natural to investigate the nonclassical symmetries associated with this equation as well. These new 
symmetry reductions also relate the forward and inverse problems as the equivalence transformations do. 
To the best of our knowledge, these type of transformations have not been studied so far for parameter 
identification problems. Since the differential structure of the equation is preserved and w is an arbitrary 
function, we shall call these transformations nonclassical equivalence transformations. While the equivalence 
transformations are found by integrating a linear PDE system, the nonclassical equivalence transformations 
are solutions of a nonlinear PDE system. Typically, a large amount of calculations are required and, hence, 
in both cases, the use of a symbolic manipulation program would represent a great advantage. Therefore, 
for this purpose, we have updated our Maple package GENDEFNC |2J (which uses the package DESDLV by 
Carminati and Vu [7]). The GENDEFNC output is the nonlinear partial differential equation system of the 



determining equations of the nonclassical equivalence transformations. Notice that GENDEFNC is based on a 
new algorithm for finding the nonclassical symmetries (see [2] for further details). Recently, Bruzon and 
Gandarias have extended this algorithm to a different case (see [H] for more information). 

As an example, we shall consider a mathematical model arising in heat conduction, namely, the nonlinear 
stationary heat conduction equation given by 

- <&v(E(x,y)Vw(x,y)) = 1 in 0,, (2) 

where the unknown function E = E(x,y) is the parameter, the arbitrary function w = w(x, y) is the data, 
and (x,y) € Q with fi C R 2 a bounded domain (here Vw = (w x ,w v ) denotes the gradient of w). The data 
function must also satisfy the Dirichlet boundary condition 

w\ dn = 0. (3) 

In 3D, the above problem is related to the heat conduction in a material occupying a domain £1 whose 
temperature is kept zero at the boundary |llj . After sufficiently long time, the temperature distribution w 
can be modeled by 

-div(E(x,y,z)Vw(x,y,z)) = f(x,y,z) in fi, (4) 

where E is the heat conductivity and / represents the heat sources. For given E and /, the forward problem is 
to find the temperature distribution w satisfying Q and (J3J). Conversely, the inverse problem is to determine 
E from ([3]) and © when w is known. While the direct problem is an elliptic PDE for w, the inverse problem 
is a linear PDE (with variable coefficients) for E. On the other hand, in the inverse problems approach, if the 
solution of the forward problem is unique for each parameter E, the parameter-to-solution map associates 
with each parameter E the forward problem solution. Since for the above problem the parameter-to-output 
map is nonlinear, ||3]) and ([3]) is a nonlinear problem from this point of view. In addition, when w — w{E), 
a new dependent variable can be introduced by using the Kirchoff transformation [251 p. 113]; if the heat 
conductivity E = 1, then the equation (J3J) becomes Poisson's equation [12;, p. 316]. For simplicity, in this 
paper, we discuss the 2D case with the heat sources / = 1. Since ([2|) can be written as 

W X E X + WyEy + E(W XX + Wyy) = ~1, (5) 

it can be easily seen that the vector field 

X = W X 8 X + Wydy - [1 + E(W XX + Wyy)] 8 E , 

generates the classical Lie symmetries associated with ([5]) considered as a PDE with w given and E unknown. 
The paper is organized as follows. The equivalence transformations related to ||SJ| are presented in 
§2 and the nonclassical equivalence transformations are discussed in §3. We show that the nonclassical 
symmetries are related to the Monge equation (j2"5| . the Monge- Ampere equation (|3ip. and the Abel ordinary 
differential equations (ODEs) of second kind (j3"8)) and (T4"6"l) . The new symmetry reductions related to ([5]) 
are given by (|3"6")) . (|55|) . l|4"4"|) . (|4T|) . and, respectively, (|4*0|) and (|4"5jl excepting the case when A(x,y) = 
(fci — k^y + k^x)j(ki + k-^x + k^y). Section 4 contains a few examples for data defined on oval and rounded 
corner domains. 

2 Equivalence transformations related to ([5]) 

Let us consider a one-parameter Lie group of transformations acting on an open set DcdxWxJ, where 
W is the space of the data functions, and £ is the space of the parameter functions, given by 

x = x + eT(x, y, w, E) + 0{e 2 ), 

y = y + eA{x,y,w, E) + C(e 2 ), 

(6) 
w = w + £$(#, y, w, E) + C(e 2 ), 

E = E + e$(x, y, w, E) + 0{e 2 ), 



where e is the group parameter. Let 

V = Y(x, y, w, E)8 X + A(x, y, w, E)d y + <P{x, y, w, E)d w + tf (ar, y, w, E)d E (7) 

be its associated general infinitesimal generator. Assume that E — E(x,y) and w — w(x,y) are both 
dependent variables in ([5]). The transformation ^ is called an equivalence transformation related to the 
PDE (J3J) if this leaves the equation invariant, i.e., w x E x + w y Ey + E(w xx + Wyy) = — 1. Note that the set of 
all analytical solutions of ([5]) will also be invariant. For this case, the criterion for infinitesimal invariance is 
given by 

pr (2) V(F)| F=0 = 0, 

with 

F (x,y,w (2) ,£ ,(1) J = w x E x +w y E y + E(w xx + w yy ) + 1, 

where pr^ 2 ) V denotes the second order prolongation of the vector field V [IH] . Notice that this prolongation is 
determined by taking into account that E and w are both dependent variables, exactly as one would proceed 
in finding the classical Lie symmetries for a PDE without arbitrary functions. The order of the prolongation 
of the vector field V is given by the highest leading derivative of the dependent variables. Applying the 
classical Lie method, we obtain the following infinitesimals 

T(x, y, w, E) = ki - k 3 y + k 4 x 

A(x,y,w,E) = k 2 + k 3 x + k 4 y ,. 

${x,y,w,E) =n(w) ( ' 

[ y(x,y,w,E)=E(2k 4 ~fi'(w)), 

where ki, i = 1,...,4 are real constants and jjl = fi(w) is an arbitrary function. Hence, the infinitesimal 
generator (0 becomes 



V = ^kiVi+V„, (9a) 



where 

Vi = d x , V 2 = d y , V 3 = -yd x + xd y , V 4 = xd x + yd y + 2Ed E , 



V M = fj,(w)d w - Efjf(w)d E - 



(9b) 



We obtain the following result 



Proposition 1 There is an infinite dimensional Lie algebra of the equivalence transformations related to 
(0) spanned by the infinitesimal generators ^. 

Therefore, the PDE |JSJ| is invariant under translations in the x-space, y-space, rotations in the (x, y)- 
space, and, respectively, scaling transformations in the (x, y, -E^-space. For instance, we find that (JSJl is 
invariant under translations in w-space if we choose /j, = const., and that the equation remains unchanged 
under scaling transformations in (w, i?)-space if we choose fJ.(w) = w. 

Furthermore, an equivalence transformation ([5]) can be used to reduce the dimension of ((SJ) by augmenting 
it with 

( T(x, y, w, E)w x + A(x, y, w, E)w y - <&{x, y, w, E) = 0, 

\ T(x, y, w, E)E X + A(ar, y, w, E)E V - $>{x, y, w, E) = 0, 

which is a first order PDE system defining the characteristics of the vector field ([7]). These relations are also 
called invariance surface conditions. 



3 Nonclassical equivalence transformations related to (|5 



Consider a one-parameter Lie group of transformations acting on an open set T> C Q x W x £, where W is 
the space of the data functions, and £ is the space of the parameter functions, given by 



x = x + s£(x, y, w, E) + 0{e 2 ), 

y = y + er){x, y, w, E) + C(e 2 ), 

w = w + ecj)(x, y, w, E) + 0(e 2 ), 

^E = E + sip(x, y, w, E) + 0{e 2 ), 

where e is the group parameter. Let the following vector field 

U = £(x, y, w, E)d x + r)(x, y, w, E)d y + <f>(x, y, w, E)d w + ip(x, y, w, E)d E 



(10) 



(11) 



be the general infinitesimal generator related to (|10|) . The transformation (|10p is called a nonclassical 
equivalence transformation (or a conditional symmetry) of the PDE ([5]) if this leaves the subset 

Sfm,4>2 = {F{x, yi w {2 \ E^) = 0,</> 1 (x,y,u> ll) , E {1) ) = 0,^(a;,y, W W, £«) = 0} 
of the set of all analytical solutions invariant, where 



4>i ■= €(%, V, w, E)w x + n(x, y, w, E)w y - cj)(x, y, w, E) = 0, 
<fc := t,{x, y, w, E)E X + n(x, y, w, E)E y - ijj(x, y, w, E) = 



(12) 



represents the characteristics of the vector field IA (or the invariant surface conditions). Here the criterion 
for infinitesimal invariance is the following 



pr 



(2 



>U{F)\p = Q > 0i=O, 02=0 = 0, 

prWW(0i)|F=o,0i=o,<fe=o = 0, 
k pr (1) W(0 2 )|F=o,0i=o,0 2 =o = 0. 

If rj ^ 0, one can assume without loss of generality that r\ = 1 (the case r\ — is not discussed in this paper), 
and, hence, (fT!?)) turns into 

w y = <p(x,y,w,E)-€(x,y,w,E)w x , ,^ 

Ey=ip(x,y,w,E)-£(x,y,w,E)E x . 

At the first step, we augment the original PDE with (I13p and eliminate all the partial derivatives of w and 
E with respect to y occurring in (|5|). Hence, by using (|13j) and its differential consequences, we obtain 



A\w xx + A 2 wl + A 3 w x E x + AiW x + A 5 E X + A 6 = 0, 
where the coefficients Ai — Ai(x, y,w, E), i = 1...6, are the following 

A 1 = E(e + i), 

A 2 = 2Etf w , 

A 3 = e + 2m E + i, 

A 4 = -& + E {& x - e y - fc w - ^e - 2^ w ) , 

A 5 = -t((f) + 2E0 E ), 

A 6 = H' + 1 - E (£cj> x - 4> y - cj)(j) w - ip4> E ) ■ 



(14a) 



(14b) 



The equation (fl4|) has been obtained by using the GENDEFNC command 



gendefnc(PDE, [w,E] , [x,y] , y, 3). 

Since A\ 7^ 0, (|14[) may be regarded as an ODE in the unknown functions w and E (with y as a parameter). 
At the second step, by using the GENDEFNC command 

gendefnc(PDE, [w,E] , [x,y] , y) 

we obtain the determining equations of the nonclassical symmetries. This is an over-determined nonlinear 
PDE system for the infinitesimals £ = £(x, y, w, E) and 4> — <f>(x, V, w, E). Among these equations, we get 

&> = 0, Ze = 0, <t> E = 0, 

which implies 

£(x,y,w,E)=A(x,y), (15) 

and 

(/>(x,y,w,E)=G(x,y,w), (16) 

where A and G are arbitrary functions. The substitution of the above functions into the remaining equations 
yields 

ip(x,y,w,E) =EF(x,y,w), (17) 

with F an arbitrary function of its arguments. By using the above relations, the determining system is 
reduced to 

F = -G W + 2 ^ A ~^\ (18) 

G x -AG y - 2AA *-Jl { f- 1) G = 0, (19) 



and 



G X X + Gyy + FyG + FGy + 2GGy W 

I 2(A x -AA y ) lnn r , Evm\ 2(AA !C +A y ) r , _„ 
~ l A 2 +l V " 1 - 7 "™ T O-j, trbj A 2 + l x ~ ' 

(A 4 - 1)A XX + AA{A 2 + l)A xv - (A A - l)A yy - 2A(A 2 - 3)A 2 X 
-4(3A 2 - l)A x A y + 2A(A 2 - 3)A 2 y = 0, 



(20) 



(21) 



where A = A(x, y), G = G(x, y, w), and F = F(x, y, w) are the unknown functions. 

Remark. The nonclassical symmetries do not leave the form of the equation |(5J) invariant. Moreover, the 
nonclassical operators (JTTJ) do not form a vector space, still less a Lie algebra, as the symmetry operators do. 
Since every classical symmetry is a nonclassical symmetry but not conversely, there exists a set of common 
solutions of the determining system of the nonclassical symmetries (nonclassical equivalence symmetries, in 
this context) and the determining system of the classical symmetries (equivalence transformations, respec- 
tively). This common solution is given by 

P(t 11 1/1 W\ — r ( x >y> w < E ) — hy-ksv+hjx 
<,yj,,y,w,J-ij A(x,y,w,E) k 2 +k 3 x+k 4 y> 

6(r 11 w E) - ®( x <y- w ' E ) - eW (921 

_ 9(x,y,w,E) _ E(2k4,-fj,'(w)) 



ip(x,y,w,E)= A \ x y y ]Z]E) 



k 2 +k 3 x+k 4 y 



where at least one of the constants ki , k^ or k^ is nonzero. 

To solve the determining equations (fl~8"l) - (f2"Tj) . we proceed as follows: the equations (fT9"|) and (|2"0)) are 
analyzed in §3.1, the equation (|2"Tj) is studied in §3.2, and the solutions of the determining equations (|18 p - (l2"l"T) 
are given in §3.3. 



3.1 The equations ([19]) and ([20]) for <j) = G{x,y,w) 
First observe that (|19p can be written in the following conservation form 

G \ f AG 



A 2 + l) \A 2 + l) °" 

' / x \ / y 

If A ^ 0, then two cases might occur: 

Case Gi. G = 0. In this case, ([TO]) and (|20p are both satisfied. 

Case G2. G ^ 0. There exists a potential function K = K(x 1 y, w) such that 

K - AG 
J^-x — A 2 +1 



<-y — A 2 + l- 

The above system yields K x = AK y whose general solution is 

K(x,y,w) =P(u,w), 
where P is an arbitrary function and u = u(x, y) is a solution of the equation 

u x = A{x,y)u y . (23) 

From these relations we obtain 

G = u y {A 2 + l)S, (24) 

where S(u,w) denotes the partial derivative P u (u,w). Thus, the solution of (fT9]) is given by ([24]) , where u 
satisfies (|23| . 

The substitution of ([24]) into ([20]) implies 

qi(S uw S - S U S W + S uu ) + q 2 S u + q 3 S = 0, (25) 

where the coefficients qi are expressed in terms of A and u y as follows 

qi =ul(A 2 + l) 3 , 

2 = Uy(A 2 + 1) [3Uyy(A 2 + l) 2 + Uy{bA X + ZA^Ay + ^A* + AAy)} , 

g.3 = u yyy (A 2 + l) 3 + u yy (A 2 + 1)[{3A 2 + 5)A X + (3A 2 + l)AA y ] 
+u y [2A(A 2 + \)A XX + (A 2 + 1)(A 2 + 3)A xy + A(A 2 + \) 2 A yy 
-2{A 2 - 3)A 2 X + 2A(A 2 - 3)A x A y + {A i - 1)^]. 

In particular, the method of separation of variables applied to (|25| implies solutions of the form S(u, w) = 
p(u)[i(w), where p = p(u) satisfies the equation 

qiPuu + q2Pu + q^p = (26) 

and \i = \i(w) is an arbitrary function of its argument. 

To summarize, the equation ([2"0]) has been reduced to ([2"5*]) . where A = A(x, y) satisfies ([2"Tj) and u = u(x, y) 
is a solution of 



3.2 The equation (J2TJ) for £ = A(x, y) 

In the following, we show that ([2"Tj) can be reduced to a Monge- Ampere equation. 

Case Al A = 

The equation (f2Tj) admits the trivial solution A = 0. 

Case AS. A = fc, with fe 7^ 

The constant function A — k, with k ^ is also a solution of (JU). 

Case A 5. A is a nonconstant function 

If A is a nonconstant function, then an alternative formulation of (j2ip is 

B(B 2 + l)fl M + {B 2 - 1){B 2 + l)B xy - B(B 2 + l)B yy 

-{W 2 - 1)B 2 - 2B{B 2 - i)B x B y + {W 2 - 1)B 2 = 

where B = (A+ 1)/(A — 1). The conservation form of the above PDE is 

B(B X + BB y )\ J B x + BB y 



(B 2 + l) 2 J x V0B 2 + l) 2 /„ 



0. (27) 



We distinguish the following two cases: 
Case A3.1. B — B(x, y) is a solution of < 
If B = B(x,y) satisfies the Monge equation 

B x + BB y = 0, (28) 

then (|2"T|) holds. Since the general solution of (j2"5|) is given implicitly by 

y-xB = v(B), 
where v is an arbitrary function, the corresponding solution of ([21 jl is 



Case A3.U. B = B(x,y) does not satisfy 

It follows from (|27|) that there exists a potential function T — T(x,y) such that 

B-r+BB,, 



T, 



__ B(B X +J 

y (s 2 +i) 



Ti B(B x +BB y ) 



(30) 



Since X^ cannot be identically zero, the above equations yield B = T y /T x and, by substituting it into the 
first equation of (j3"0|) we have 

T x T y T xx - (T 2 - T 2 ) T xy - T x T y T yy + {T 2 + T 2 f = 0. 
By using the following Legendre transformation [251 p. 353] 

7i(a, b) + T(x, y) — xa + yb, 



where 



where T x = a, x = 7i a , T y = b, and y = TL^, the above PDE turns into the following Monge- Ampere equation 

[a z + b z ) (a z + b z ) [a z + b l ) 

Furthermore, this can be reduced to the Monge- Ampere equation 

VaaVn - Vl = — l (31) 

(a z + b z ) 

V(a, b) = H(a, b) - - arctan f|j . (32) 

Exact solutions for particular Monge- Ampere equations have been extensively analyzed in 22J. The Monge- 
Ampere equation (|3"Tj) may be included in Case 17 [2H p. 458] or Case 20 [22, p. 460]. 

3.3 Solutions of the determining equations ( 1181) — ( 12TT) 

We distinguish the following four cases: 
Case 1. A = and G = 

Proposition 2 If A = and G = 0, the infinitesimal generator ill]) becomes hi = d y , which implies the 
invariance of the equation (5p with respect to translations in the y-space. 

This nonclassical symmetry reduction is an equivalence transformation that can be obtained from (|22p 
for ji = 0, &2 = 1, and fej = 0, where % = 1, 3, 4. 

Case 2. A = and G =£ 

The PDE HU) takes the form G x = and it follows that G = H(y,w). After substituting it into (JTSJ) and 
(|2H)l . we get F = — ff„, and, respectively, 

Hyy + HHy W ~ HyH W = 0, (33) 

where iJ is an arbitrary function. Notice that H ^ in the above equation. We distinguish the following 
cases: 

Case 2.1. H = /J,(w), where /i is an arbitrary function 

Indeed, H = /J,(w) is a particular solution of (|3"3)) . 

Proposition 3 If A = and G = fx{w), where fi is an arbitrary function, the infinitesimal generator ill]) 
turns into 

U = d y + f i(w)d w ~ f i'{w)Ed E . (34) 

The above nonclassical equivalence transformation is in fact an equivalence transformation and corre- 
sponds to the case &2 — 1) an d k, = 0, where i = 1, 3, 4 in 



Gase 2.2. H y =£ 

The PDE |33|) can be written in the following conservation form 



After introducing the potential function g = g(y,w), we get 

{ 9w = —jj^. 

Clearly, g w ^ 0. After eliminating H y in the above system, we obtain H = —g y /g w - Next, substituting it 
into the second equation, the following PDE results 

9y(9y ~ v)w - 9w(g y - y) v = 0. (35) 

The following two cases occur: 

Case 2.2.1. g y = y 

It results g(y,w) = y 2 /2 + h(w), where h is an arbitrary nonconstant function (otherwise, g w = 0). Since 
H = —g y /g w , we get H — yfi(w), where fi = —1/h'. 

Proposition 4 If A = and G = y/i(w), where fj, is an arbitrary function, then the infinitesimal generator 
Ul\) becomes 

U = d v + yn{w)d w -yn'{w)Ed E , (36) 

The nonclassical symmetry generated by (|36[) is a new symmetry reduction for ([5]) which cannot be 
obtained from 



Case 2.2.2. g y - y =£ 

Since the equation (|35|) is the Jacobian of the functions g and g y — y, there exists a function a such that 

g y = y + a(g). (37) 

In the above relation, w is viewed as a parameter. The equation (|37[) can be written as the following Abel 
ODE of second kind 

(y + a(g))^- = l, 
dg 

for y — y(g). The canonical substitutions z = a{g) and v = y + Q.(g) reduce the above ODE to its canonical 
form 

vv'-v^Piz), (38) 

where v — v(z) and j3 = l/(a' o a~ l ). A collection of the known cases of solvable Abel ODEs of the form 
(|46)l is presented in [21j pp. 107-120] and new results can be found, for example, in [8]. Each of these ODEs 
corresponds to a nonclassical equivalence transformation related to ([5]) . 



Proposition 5 // A = and G = —g y /g w , where g is a solution of Ij37\ ), then the infinitesimal generator 
ill]) turns into 

U = d y - g y /g w d w + E(g y /g w ) w d E , (39) 

Since g satisfies ([35|) . the above vector field generates new symmetry reductions related to ([5]) that are 
not equivalence transformations. 

Case 3. A =£ and G = 

By HID, we obtain F = 2(A X - AA y )/(A 2 + 1). 

Proposition 6 If A ^ is a solution of the equation H21\) and G = 0, then the nonclassical infinitesimal 
generator ill]) becomes 

U = A(x, y)d x +d y + 2(Ax A ~^ v) Ed E . (40) 
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Case 3.1. A — k, where k ^ 

Proposition 7 If A = k, where k ^ is a constant, the infinitesimal generator |^0[ ) rewrites as hi = 
kd x + d y . 

Replacing (i = 0, fe = 1, and ki — with i = 1, 3, 4 in (|22[) . we obtain the symmetry reduction generated 
by the above vector field, and, therefore, this is an equivalence transformation. 

Case 3.2. A is a nonconstant function 

If A is a nonconstant function, then we distinguish two subcases: 

Case 3.2.1. A is given implicitly by 1^5]). 

Case 3.2.2. The equation for A is reduced to the Monge- Ampere equation pTj) . 

In the above two cases, the nonclassical operator (|4"0)) generates new symmetry reductions for ([5]) excepting 
the case when A(x, y) — {k\ — k^y + k^x)j(ki + k$x + k$y). 

Case I A =£ and G ^ 

Case 4-1- A= k, where k ^ 

Without loss of generality, we consider the particular solution u(x.y) = kx + y of the equation (|23|) . Since 
My = 1, the relation (124[) yields G = (fc 2 + l)^, where S is a nontrivial solution of the equation (f2~5|) which 
rewrites as 

iutl + bb uw — O u D w = 0. (41) 

The following cases might occur: 

Case 4.1.1. S = /J,(w), where [i is an arbitrary function 

In this case, ([4"Tj) is satisfied and (fT5| implies F = —/j/(w). 

Proposition 8 If A = k with k ^ and G = (A; 2 + l)//(w), where /1 is an arbitrary function, then the 
nonclassical infinitesimal generator Ml]) becomes 

U = kd x + d y + (k 2 + l)/j,(w)d w - (k 2 + l) fJ /(w)Ed E . (42) 

The above nonclassical operator generates an equivalence transformation that can be obtained from (|22|) 

for fci = k/{k 2 + 1), fc 2 = l/(k 2 + 1), fc 3 = 0, and fc 4 = 0. 

Case 4.1.2. S u ^ 

Since the PDE (|4Tj) can be written in the conservation form 

there exists a potential function Q = Q(u,w) such that 

{^ u = ST' 

Eliminating S in the above system, we obtain 5* = —Qu/Qw, where 

Qu(Qu ~u) w - Q W (Q U - u) u = 0. (43) 

We distinguish the following two cases: 
Case 4-1-%-a- Qu = u 

In this case, Q(u, w) = u 2 /2 + p(w), where p is an arbitrary nonconstant function (otherwise Q w = 0). With 
the aid of S = —Q u /Qw, we get S(u, w) = u^jl{w), where // = —1/p' and u(x, y) = kx + y. 
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Proposition 9 For A(x,y) — k (k ^ 0) and G(x,y,w) = (k 2 + l){kx + y)ji{w), where /i is an arbitrary 
nonconstant function, the nonclassical infinitesimal generator \11\) rewrites as follows 

U = kd x + d v + {k 2 + \){kx + y)fi(w)d w - (k 2 + \){kx + y)n'{w)Ed E . (44) 

The nonclassical operator (|44[) generates new transformations that are not equivalence transformations 
related to ([5]). 

Case 4.1. 2.b. Q u - u ^ 

Observe that (|4"5|) is the Jacobian of Q and Q u — u. Thus, there exists a function 7 such that 

Q„ = u + 7(<9)- (45) 

In the above ODE, w is viewed as a parameter. Similar to Case 2.2.2, the equation ([45]) can be reduced to 

/ ,^, du 

(«+ 7(g» 55 = 1, 

which is an Abel ODE of second kind for u — u(Q). Moreover, after using the substitutions s = j(Q) and 
V = u + j(Q), the above ODE can be reduced to the canonical form 

W - V = 6(s), (46) 

where 6 = 1/(7' o 7 _1 ). This is an Abel equation of second kind for the unknown function V = V(s). For 
each solution of (|46|) . we obtain a nonclassical symmetry for ([5]). The solvable Abel ODEs of the form (|46]) 
that are known so far are listed in [23 pp. 107-120]. More recent results can be found, for instance, in [S]. 

Proposition 10 If A[x,y) = k with k 7^ and G(x,y,w) = — (k 2 + l)Q u /Qw> where u(x,y) — kx + y, and 
Q = Q(u,w) satisfies J^ffi ), then the nonclassical infinitesimal generator ill]) is given by 

U = kd x + d y - (k 2 + 1)Q U /Q w d w + (k 2 + l)E(Q u /Q w ) w d E . (47) 



Since Q satisfies (|43|). the above vector field generates new symmetry reductions related to (|5|) that are 
not equivalence transformations. 

Case 4-. 2. Suppose A is a nonconstant function and G ^ 0. 



Proposition 11 If A is a nonconstant function satisfying the equation A21\) and G = u y (A 2 + 1)S, where u 
satisfies (23\) and S is a nonzero solution of i25\), the nonclassical infinitesimal generator \11\) is written as 

U = Ad x + d y + u y (A 2 + l)Sd w + E (-u v (A 2 + 1)S W + 2( ^ 2 ~ + A 1 Ay) ) d E . (48) 

Two cases might occur: 

Case 4-2.1. A is given implicitly by (129]) . 

Case 4 ■"2. 2. The equation for A is reduced to the Monge- Ampere equation (f3Tj) . 

In the above two cases, the vector (|48p generates new symmetry reductions that cannot be obtained from 



excepting the case when A(x, y) = (fci — k$y + kix)/(k2 + k^x + k^y). 
We thus have found new symmetry reductions related to ([5]) and have shown that these are given by (j36[) , 
(|39| . (|44]1 . (|47]) . and, respectively, (j40| and (JUJ excepting the case when A(x,y) = (fci - k 3 y + k±x)/(k 2 + 
k^x + k 4 y). 
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4 Examples 

The nonclassical symmetries related to ([5]) yield classes of data and suitable domains for which the dimension 
of the problem can be reduced. Since the homogeneous Dirichlct boundary condition ([3]) is imposed, the 
data w must satisfy (I13[) and the boundary dQ must be " compatible" with the symmetry reduction as well. 
To illustrate this approach, we discuss a few examples of data modeled by the functions 



w(x, y) = —y 2 + W(x) 



(49) 



which are invariant with respect to the nonclassical infinitesimal generator 

U = d y - 2yd w , 

which is obtained from (3.27) by setting (J,(w) — —2. Substituting ([50)1 in ([l3"|l. we get w y - 
This implies E(x,y) — Eq(x) (i.e., the heat conductivity is constant along the lines x = 
given by (|4D|) . Therefore, we can consider circles, ellipses, generalized Lame curves x 2p 



(50) 

— 2y and E y = 0. 

const.) and w is 

-y 2 = 1 (p>2), 



Granville's egg curve y 2 x 2 = (x — b)(l — x), where 6 ^ 0,1, or elliptical curves - in particular, Newton's 
egg curve y 2 = (x 2 — l)(x — a), where a^±l. In this case, as it is shown below, the parameter cannot be 
determined at the points (xq, y) for which W'(xq) — 0. Despite this fact, additional information about the 
parameter on the specified domain can be obtained. 

Example 1. (Newton's egg curve) Suppose dfl = {(x,y) : y 2 = (x 2 — l)(x — 3), x <G [— 1, 1]} (see Figure 
[1|). For the data 

w(x,y) = -y 2 + (x 2 - l)(x - 3), 

the equation ([5]) is reduced to the ODE (3x 2 — 6x— l)E' (x) + (6x — 8)Eq(x) — — 1, whose general solution is 



C-E 1 (x) 



.^,^-d 



for x ^= 1 — 2/-\/3, where E[(x) = exp [y/3/3 arctanh ((.t — l)v3/2)l . Figure O shows the parameters satis- 
fying the conditions E(— 1) = 0.2, and, respectively, E(l) = 0.2. The parameter E cannot be determined for 
.T = 1 - 2/V3. 





Figure 1: The boundary 9fi = {(x,y) : y 2 — (x 2 — l)(x — 3) = 0, x G [—1,1]} and the data w(x,y) = 
—y 2 + (x 2 — l)(x — 3) discussed in Example 1. 
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-0.8 -0.6 -0.4 -0.2 



Figure 2: The parameter E(x,y) = Eq(x) corresponding to data discussed in Example 1. 

Example 2. (Granville's egg curve) Consider dVt = {(x, y) : x 2 y 2 = (x — 3)(1 — x), x G [1, 3]} (see Figure 
[3]) and the data given by 

2 , (z-3)(l-x) 



w{x,y) 



-y 



The reduced ODE 2x(2a; - 3)E' {x) + 2(x i - Ax + 9)J5q(x) = a; 4 implies 



£(x, y) = £ (z) = L (2a , _ 3) V; 16 exp 



— (4a; 2 + 9a; + 27) 
24 



which is defined for a; + 1 1.5, where E[(x) = 0.5 (2a; - 3) 27/16 exp [|j(4a; 2 + 9a; + 27)] . Similarly, as in 
Example 1, the parameter cannot be determined for x = 1.5. See Figure H for E(l) — 0.2 and, respectively, 
E(3) = 0.6. 





Figure 3: The boundary dfl — {(x,y) : x 2 y 2 = (x— 3)(1— x), x e [1,3]} and the data w(x,y) = y 
considered in Example 2. 



_ „, 2 (x-3)(l-x) 
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Figure 4: The parameter E(x,y) = Eq(x) corresponding to data considered in Example 2. 

Example 3. (Generalized Lame curve) For dil — {(x, y) : x A + y 2 = 1, x e [— 1, 1]} (see Figure[5]) and the 
following data 

w(x,y) = 1-x 4 ~y 2 , 

the PDE © becomes 4x 3 E' (x) + 2(6x 2 + !) = !. It results 



*(*) = i 



—^[ 


' n 


+^1 


8x 3 ( 


K 2xj 


a; 3 J 



exp 



1 
4x 2 



that is defined for x ^ (here erf denote is the error function). Therefore, the parameter cannot be estimated 
at x = 0. See Figure [B] for the cases Eq(— 1) = 0.25, and, respectively, -Eq(I) = 0.25. 






Figure 5: The boundary dil — {(x,y) : x 4 + y 2 = 1, x e [— 1, 1]} and the data w(x,y) = l — x A — y 2 discussed 
in Example 3. 



Example 4- (Rounded corner rectangular domain) Suppose 90 = {(x, y) : y 2 = — 8x 4 +x 2 +2, x € [—a, a]}, 
where a = V 1 + V65/4 (see Figured]). For the data 

w(x,y) = 2 + x 2 ~8x i -y 2 , 
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Figure 6: The parameter E(x,y) = Eq(x) corresponding to data discussed in Example 3. 





Figure 7: The boundary 9f2 = {(x,y) : y 2 = —8a; 4 + x 2 + 2, x e [—a, a]}, where a = \J\ + \/65/4 and the 
data w(x, y) = 2 + x 2 — 8a; 4 — y 2 considered in Example 4. 



we obtain the following reduced equation 2a;(16a; 2 — l)E' Q (x) + 96x 2 Eq(x) — 1. Thus, 



E (x) = 
for x S [-a, -0.25) U (0.25, a] , and 

E (x) = 



C x + \/16.x 2 - 1 + arctan (Vl6a; 2 - l) 
2(16a; 2 - l) 3 / 2 



-1/2 



C 2 + Vl - 16a; 2 + arctan (Vl - 16a 
2(1 - 16a; 2 ) 3 / 2 



-1/2 



for a; £ (—0.25,0.25). E 1 cannot be determined for x = ±0.25. Figure [8] shows the graph of the parameter 
for d = C* 2 = 0. 



1G 




Figure 8: The parameter E(x, y) = Eq(x) corresponding to the data discussed in Example 4. 

5 Conclusion 

In this paper we point out another systematic way of finding classes of symmetry reductions related to 
parameter identification problems of the form |T]). Similar to the study in [3J, we emphasize that the 
geometrical significance of the nonlinearity occurring between the data and parameter in (fTJ) can be reflected 
by the group analysis tool. Seeking different shapes for domains on which the dimension of the problem can be 
reduced is not an easy task. Therefore, in this paper, we discuss the nonclassical equivalence transformations 
related to ([5]). Briefly, to determine these transformations, the data is considered as a dependent variable 
as well as the parameter, and the nonclassical method (due to Bluman and Cole) is applied to the studied 
equation. In general, for a known data w, one should check the invariance of this function in (|13p . where £, 
<j), and tp are discussed in §3. Next, from the second equation in (|13p . the form of the parameter should be 
obtained in terms of the invariants of the symmetry reduction. At the end, substituting E and w into iJSJ , 
the dimension of the model should be reduced by one. 
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